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ABSTRACT The occurrence of scissors modes in crystals that have deformed ions in their cells 
has been predicted some time ago. The theoretical value of their energy is rather uncertain, however, 
ranging between 10 and a few tenths of eV, with the corresponding widths of 10“^, 10“® eV. Their 
observation by resonance fluorescence experiments therefore requires a photon spectrometer covering 
a wide energy range with a very high resolving power. We propose and discuss a new experiment 
in which such difficulties are overcome by measuring with a SQUID the variation of the magnetic 
field associated with the excitation of scissors modes. 

PACS numbers: 75.10.-bGeneral theory and models of magnetic ordering, 71.10.-wTheories and models of 
many-electron systems 


I. INTRODUCTION 

In a series of papers it has been suggested that deformed atoms in crystal cells might have collective excitations 
called scissors modes nHa. These are states in which two particle systems move with respect to each other conserving 
their shape. They were first predicted to occur in deformed atomic nuclei [1] by a semiclassical Two-Rotor-Model 
(TRM) in which protons and neutrons were assumed to form two interacting rotors to be identified with the blades of 
scissors. Their relative motion (Fig.l) generates a magnetic dipole moment whose coupling with the electromagnetic 
field provides their signature. 

After their discovery [S] in a rare earth nucleus, and their systematic experimental and theoretical investiga¬ 

tion [6] in all deformed atomic nuclei, scissors modes were predicted to occur in several other systems including metal 
clusters [7], quantum dots [S], Bose-Einstein [5] and Fermi m condensates and crystals m- In all these systems 
one of the blades of the scissors must be identified with a moving cloud of particles (electrons in metal clusters and 
quantum dots, atoms in Bose-Einstein and Fermi condensates, individual atoms in crystal cells) and the other one 
with a structure at rest (the trap in Bose-Einstein and Fermi condensates, the lattice in metal clusters, quantum dots 
and crystals). The dynamics of such systems can then be described by a One-Rotor-Model (ORM). 

More recently the TRM has been used m also for single domain free nanoparticles m- These objects consist of a 
magnetic structure, called macrospin, that rotates with respect to a non magnetic lattice. They have been represented 
as a couple of rigid rotors, one associated with the nonmagnetic lattice, and the other one, with a spin attached, with 
the macrospin. If the nanoparticles are not free, but stuck in rigid matrices, the dynamics is determined by a ORM 
with spin, representing the macrospin. 



FIG. 1: Scissors Mode in atomic nuclei: the proton and neutron symmetry axes precess around their bisector 
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FIG. 2: In a) an atom is in its ground state in a cell while a photon is incoming with circular polarization and momentum 
parallel to the axis of the cell. In b) the photon has been absorbed transferring its angular momentum to the atom, which 
processes around the symmetry axis. A photon with momentum orthogonal to the cell axis cannot be absorbed, because the 
atom cannot absorb its angular momentum since it cannot rotate around any axis orthogonal to the cell axis. 


Scissors modes have been observed not only in atomic nuclei, but also in Bose-Einstein condensates m- Moreover 
nanoparticles stuck in a rigid matrix were studied in detail by a ORM, and the magnetic susceptibility was found 
compatible with a vast body of experimental data and in some cases the agreement was surprisingly good m- So 
we can fairly state that the Two-Rotors or the One-Rotor models are relevant to disparate physical systems involving 
energies and sizes at very different scales. 

Scissors modes in crystals have also been studied in the framework of a ORM in which an ion is regarded as a 
rigid body, which can rotate around the axes of its cell under the electrostatic force generated by the ligands. We 
considered crystals with uniaxial [T] and cubic symmetry [5]. In the hrst case the precessing ion was treated as one 
ellipsoidal rotor, in the second case as the body obtained by superimposing three ellipsoids at right angles. In the 
presence of uniaxial symmetry the photoabsorption cross section is characterized by a linear dichroism [1] (Fig.2) due 
to the fact the ion can rotate around the symmetry axis of the cell, but not around the other ones. 

In the first paper on crystals [T] it was proposed to search for scissors modes with a resonance fluorescence experi¬ 
ment. In such an experiment the dichroism can help to eliminate the background. 

Then it was considered the case of ions characterized by Spin-Orbit Locking, namely ions in which the spin-orbit 
force is so strong to lock the total spin to the density profile in such a way that they rotate together. For such a case 
a different experiment was proposed, inspired by the experiments on Bose-Einstein condensates m- In the latter 
ones one induces oscillations of the condensate giving a sudden twist to the magnetic trap. In the case of ions with 
spin-orbit locking one induces oscillations of the ions in the cells by an impulsive change of an external magnetic field 
and observes the photons emitted when the ion goes to the ground state. 

The estimate of the energy of scissors modes is affected by a large uncertainty, ranging from 10 to a few tens of 
eV. This makes the experiments difficult, because one have a spectrometer covering a large x-ray energy range with 
a resolving power up to 10^. Such a high value is due to the long lifetime of the scissors mode and could explain why, 
as far as we know, no hint of scissors modes shoved so far in the experimental measurements of spectra in crystals. 

We therefore propose to excite scissors modes by an intense source of x-rays, and instead of detecting the photons 
emitted in the decay, to measure the variation of the magnetic field associated with such excitation. To this purpose 
we can use a Superconducting Quantum Interference Device (SQUID), that is suitable to measure extremely small 
magnetic helds. 

The restoring force acting on the ions in our model is mainly due to the interaction of their electric quadrupole 
moment with the electric quadrupole moments of the ligands. Scissors modes will then be realized easier in compounds 
in which there are highly deformed ions, and natural candidates are compounds containing Rare Earths (RE). Moreover 
the /-electrons of the Rare Earths are compact and for this reason they can be reasonably assimilated to a rigid rotor. 

We will consider crystals and gels in which Rare Earths are dispersed. The second case is very interesting, because 
the RE ions are well separated, so as to make completely justified our basic approximation that they are non interacting 
with one another. It will turn out, unfortunately, that our type of experiment with gels looks below the limit of present 
possibilities. 

We end this Introduction by noticing the many reasons of interest of Scissors Modes. On general grounds for every 
system they should affect the dispersive effects in the channels with their quantum numbers, which are J'^ = I^. 
For given systems they provide specihc pieces of information. In nuclear physics they are related to the superfluidity 
of deformed nuclei, in Bose-Einstein condensates provide a signature of superfluidity, in metal clusters they are 
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predicted to be responsible for paramagnetism, in nanoparticles they give information about the moment of inertia 
of the macrospin. The knowledge of their properties should be of some consequence also for crystals, because the 
magnetic anisotropy of these systems is at the origin of many interesting technological applications including magnetic 
storage devices and sensors, spin-torque nano-oscillators for high-speed spintronics and spin-optics M- In particular 
one could establish or rule out the existence of Spin-Orbit-Locking in RE atoms. 

We will start by describing, in Section II, the ORM that we will use to describe the dynamics of a deformed atom 
in a crystal cell or in a gel. Next we will discuss in Section III the estimate of the parameters entering our model. 
Finally in Section IV we will evaluate the magnetic field generated by photon absorption to be measured by a SQUID 
and will present our conclusions in Section V. 


II. THE ONE-ROTOR-MODEL 


In the Single-Ion Model of magnetic anisotropy each rare earth ion is assumed to be independent from the others. In 
the study of the magnetic properties only the motion of the 4/ electron system as a whole is considered, disregarding 
its excitations. In other words this system is treated as a rigid rotor with spin. We introduce the principal frame of 
inertia of the ion with axes 77 , ^ and study its dynamics with respect to a frame of reference fixed with the cell with 
X, y, z-axes parallel to the cell axes. We restrict ourselves to a rotor with axial symmetry and assume the symmetry 
axis parallel to the C-axis. In quantum mechanics rotations around the symmetry axis are forbidden, and therefore 
the wave functions must satisfy the constraint 

C-LV'=0 (1) 


where L is the angular momentum the ion. Its components in the cell frame become then those of a point particle. 
Assuming the symmetry axis of the ion along the Qaxis the kinetic energy operator is 


T = 





— cot 9 


de 



( 2 ) 


where 9 is the angle between the z,('-axes, = —ih, is the z-component of the orbital angular momentum, and I 
the moment of inertia with respect to the and Ty-axes. For the potential we assume 


V = -Csm^9 
2 


( 3 ) 


where C is a restoring force constant. 

The scalar product for such a system remains defined by the measure of integration over the Euler angles, which 
for functions independent of 7 is 


P‘2'K nTT 

(V^il^/)=/ dip d9 sin 9 Ip* {9, p)ijjf {9, p) = 6if (4) 

Jo Jo 

The exact eigenfunctions and eigenvalues are known m, but in the cases of physical interest the quadratic approxi¬ 
mation is sufficient 


T 

V 




lf2C9^ 

I/2C(7r- 


9 d9^ 9^^^ )' 


0 < 6» < TT 


, 0 < 0 < 7 r/ 2 , 

9 , 7r/2 < 9 < TT. 


( 5 ) 


We recognize that this Hamiltonian, in each of the 0-regions above, is exactly that of a two-dimensional harmonic 
oscillator, provided we identify 0 with the polar radius. As shown by the following estimates, the fall off of the wave 
function is so fast that we can extend without any appreciable error the integral over 0 up to infinity. 

In the quadratic approximation, in the region 0 < 0 < 7 r /2 the eigenfunctions are 

'0n,iJ0,<p) = e*'"‘^x„,|j^|(0) (6) 

where 


( 7 ) 




4 


L 




are Laguerre polynomials and 


The eigenvalues are 




1 2n\ a2 _ ^ 


OoV {n+\h\)V V^' 


En^i^ =huj{2n+ \lz\ + 1), (n = 0, 1, 2 ...) 


( 8 ) 

( 9 ) 


where 


( 10 ) 

The first excited states have n = 0, Iz = ±1, and as we will see, they are the only states strongly coupled to the 
ground state by electromagnetic radiation. As shown in Fig.2, they describe the precession of the atom at an angle 
0 ^ 00 around the axis of the cell. 

If the ion is spinless its wave function is invariant under the inversion of the orientation of its symmetry axis, that 
is not observable. A discussion of the consequences of this symmetry can be found in [16], but they will be ignored 
here because most RE ions do have a nonvanishing spin. 

The dynamics is determined by the magnitude of the zero point fluctuations 9q. When 6q ^ 0 the axis of the 
rotor can be assumed to lie along the z-axis, the direction of easy magnetization, and its zero-point fluctuations can 
be ignored. For 0o ~ 1 the zero-point fluctuations cannot be neglected, but the rotor is still polarized within an angle 
of order 0^. For 6 q » 1 there is no polarization at all. 

We quote the values of Oq in some systems: in the atomic nuclei [I] of the rare earths, 0g « 10“^; in the crystal 
LaMnOa, 0.1 < 0g < 0.5, depending on how the moment of inertia and restoring force constant are evaluated. Note 
that in the mentioned compound the deformed ion is Mn. 



III. ESTIMATES OF THE PARAMETERS 

For numerical evaluations we need estimates of the moment of inertia I and of the restoring force constant C of 
the RE ions. 

For the moment of inertia we should assume, consistently with our One-Rotor Model, the expression appropriate 
to a rigid body 


I=‘^meZ{rl) ( 11 ) 

where is the electron mass, (r^) the mean square radius of the electrons that contribute to the moment of inertia 
and Z their number. We must then establish which are the electrons that effectively take part in the rotation. Indeed 
one might include or exclude in the evaluation of I those electrons whose distribution has a spherical shape |3] . For 
a RE this means to include or exclude all the electrons of the inner shells and the 4/-electrons that have a spherical 
density. In the presence of Locking, however, the 4/-electrons that have a spherical distribution have a non-vanishing 
spin and therefore rotate with the electrons which determine the charge deformation, so that the ambiguity is restricted 
to the inner shells. If we assume that only the 4/-electrons participate in the rotation, we can set (r|) = (rlf) « ag 
(ag « 0.5 A) and Z - Zif « 10 if we consider, for instance, Dy, Ho and Er, getting 

-«8eE. (12) 

The main contribution to the restoring force comes from the interaction between the electric quadrupole moment of the 
RE ion and the electric quadrupole moment of the ligands. So the expression of C depends on the specific structure. 
One example of its calculation can be found in m, but here we do not want to restrict ourselves to particular cases. 
Therefore we will avoid the evaluation of C and parametrize the energy of the scissors modes according to 



IDS 


-2 

0 


eE. 


Es = h 


(13) 
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We will assume that 9q must be smaller than 1 if the RE ion must be polarized, and in specific cases we indeed found 
el « 0.5. 


It is also interesting to have an estimate of the resolution power for our experiment. The expression of the width 
of the scissors mode was determined in Refs. nia 


r « -a — 


1 


(WeC^) 


2\2 0 


0 , 7 ® > 1.6 • 10-^°eV. 


(14) 


Notice the strong dependence on 9q. With our parametrization of the energy we get 


Es 

Y' 


4 

3a 


ITleC^I 


06 = 6 • lO®0g < 7 • 10^ 


(15) 


We must notice that the above expression is not entirely consistent, because while the expression of the scissors mode 
energy has a general validity, the width was evaluated for a specific structure of the a crystal cell. Nevertheless we 
think it should be sufhcient to give an order of magnitude of the resolving power. 


IV. ESTIMATE OF THE MAGNETIC FIELD GENERATED BY PHOTON ABSORPTION 


We assume that all the crystal cells have their easy axis parallel to the z-axis of the laboratory frame. This situation 
is realized in single crystals and can be easily achieved by applying a convenient magnetic field. Otherwise an average 
over the cell axes should be performed as in m- We also assume that the photons are circularly polarized (left or 
right is not relevant). The change of the magnetic moment of the target is given by the probability to excite a scissors 
mode times the scissors mode magnetic moment minus the ground state magnetic moment. It can be seen that the 
spin contribution is 0g times smaller than the orbital one and therefore it will be neglected. The orbital magnetic 
moment in the ground state vanishes and then 

A = ui tfiux Na < tpoh I > e-B (16) 

where tfi^x is the duration of the photon flux, Na the effective number of atoms in the sample and the Bohr 
magneton and w is the photon absorption probability per unit time per unit frequency, that is given by 

= 47r^aV(w)^^^|(Me(w))/ip . (17) 

In the above equation rrie is the electron mass, a the fine structure constant, c the velocity of light, N{uj) the number 
of photons of energy lo per unit volume and M(oj) the magnetic dipole moment matrix element 

(Y)f, =-^ < > (18) 

where 

and e = ±1 is the photon polarization. 

Using the expression of the wave functions, for the quantum numbers f = {0,m),i = (00) 


(M, 


lOJ 


<^)0h,00 


2V2C0O 


( 20 ) 


so that 


'4^oiz\ ^— A? A — ii 


= —a 


^phot 


2 02 " Aw S',. 


( 21 ) 


( 22 ) 
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where Ag is the Compton wave length of the electron, u) the photon frequency, Aw the frequency spread, Nphot the 
number of incoming photons per unit time and Sgource the cross section of the photon beam. 

The effective number of atoms is given by 


— P ^phot Ssamp 


(23) 


where p is the number of atoms per unit volume, Ssamp the sample surface and Xphot the photon rms path. The above 
formula holds provided Xphot is smaller than the sample thickness. 

We assume the bandwidth proportional to the photon frequency according to 


Aw — f lOphot — f ^sciss 


(24) 


so that finally 


I A 


1 1 


^ 02 r Xlphot if lux Q Ag Xphot P P^B • 

^ t/Q J '^source 


A reasonable assumption is that 


5.. 


■^samp 1 


/= 10 


-3 


SO that 


I A 



o2 i^phot i Ag Xphot P 10 PB ■ 
^0 


(25) 


(26) 


(27) 


The strength of the magnetic field along the z-axis generated by this variation of the magnetic moment at a distance 
r is 


B, 


^Pz 

— 


where /xq is the vacuum magnetic permeability. 
Next we assume 


Nphot « 10“sec ^ 


(28) 


(29) 


el < 0.5. (30) 

To proceed further with numerical estimates we must distinguish the case of a dense RE compound sample as a crystal 
from the case of gels containing highly diluted RE atoms. 

All the following numerical estimates will be done in S.I. units. Then 

Aio = 1-25 • 10-®iJ/m, /XB = 9.3-10-24 J/T (31) 


and 


Ag = 4 • 10-43 to . (32) 

A. Dense matter 

We assume one RE atom per cell and a cell volume of (3.5A)^ = 4.3 • m®, so that pnE = 2.3 • 10®® to”®. 

We then have 


I A /x^l — 7 • 10 tfiux Ag Xphot P Pb 


(33) 


and 


B — 7 • 10 tfiux ^phot P'^ PoPB • 


( 34 ) 
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Setting 


Ap/iot ~10 "‘to, r~10 ^TO 


(35) 


we get 


S~4 - 10 in seconds. (36) 

Assuming an acquisition of one hour we should measure a magnetic moment of the order of 10“^ T, which is lower 
then the Earth’s magnetic field on the surface whose magnitude ranges from 20 to 100 /xT. However we can increase 
the signal reducing the distance r of the detector and increasing the acquisition time and the flux intensity. 

Examples of suitable compounds for the proposed experiments are R 2 Eei 4 B and R 2 Eei 7 N 3 in which the RE ions 
can be Dy, Ho and Er. We remind that the spin contribution to the variation of the magnetic moment is 6q times 
smaller than the orbital contribution. Since all the quoted RE’s have Z 4 / « 10 and gJ « 10, if Oq is not sufficiently 
small the spin contribution cannot be neglected. 

Einally we remark that the ions Gd^’*', La^+, being spherical, do not have an orbital response. They can then be 
used to exclude, by comparison, spurious effects. 


B. Soft matter 


The difference with respect to a crystal is due to the difference of RE density and photon rms path 

{PRE^gel {,^phot)gel 

^gel — -^crystal 7 ^ TT T 

\PREjcrystal \^photjcrystal 

If X is the percentage in weight of RE atoms in a gel made of Ngei atoms of mass vrigei 


(37) 


(PRE)gel ITlRE = X PgslRlgel ■ 

Assuming mgei/niRE ~ 0.1, pgei ~ 0.5 • 10^^ m~^ 


(38) 


{PRE)gel « 3^ ' 0.5 • 10 ® {PRe) crystal ' 


(39) 


With a percentage in weight of x=0.05 we will have a magnetic moment B « 10“^^ tfiux T . Considering that 10“^^ 
T is the minimum value of the magnetic field actually detectable by existing SQUID devices we should imagine an 
acquisition of tens or even hundreds of hours or reduce significantly the distance and increase the flux intensity to be 
able to detect a signal from such a diluted RE system. 


V. CONCLUSION 

The previous experiments designed to observe scissors modes in dense matter are affected by two orders of difficulty. 
The first one is the detection of photons in a wide energy spectrum, due to the uncertainty in the theoretical value of 
the energy of these excitations, that ranges between 10 and a few tens of eV. The second one is the need of a resolving 
power of the order of 10^, much larger than that routinely used, that is of the order of 10^ — 10^. These difficulties 
might also explain why, if scissors modes exist, they have not been noticed until now. 

We propose here an experiment in which, instead of detecting the decay photons, one has to measure with a sensitive 
SQUID the magnetic field associated with the excitation of scissors modes induced by photon absorption. For such 
an experiment we only need a sufficient flux of photons and a spectrum of the source sufficiently wide to cover the 
theoretical range of the scissors modes energy. 

There are some features of the present experiment that should help to eliminate spurious signals 

1 ) there should be no magnetic response with a source of linearly polarized photons 

2 ) the intensity of the variation of the magnetic field should decrease as the inverse of the cubic distance between 
the SQUID sensor and the sample 

3) there should be no magnetic response from spherical RE ions. 

4) because our calculations are performed assuming that the cells axes are aligned, the induced variation of magnetic 
field should be proportional to the square of the cosine of the angle between the direction of the photons and the easy 
axes of the cells. 



Alignment of the axes is realized in single crystals, and can also be achieved by applying an external magnetic field. 
As far as RE dispersed in gels matrices are concerned, it appears possible to synthesize samples with a signihcant 
alignment. 

If the direction of the cells axes are disordered, we must perform an average as done in m for nanoparticles. 
Obviously such an average does not change the order of magnitude of the results. 

We must distinguish now between dense compounds and light systems as RE ions dispersed in a gel matrix. 

In the first case with the available photon fluxes in a synchrotron radiation facility in the ultraviolet energy range, 
we estimated the variation of the magnetic field to be of the order of 10“^ Tesla. This is to be compared with the 
typical precision of measurement of the terrestrial magnetic field that is of the order of one nanoTesla. 

In the second case, with a small amount of RE ions dispersed in a gel matrix, our estimate yields a magnetic 
field whose magnitude is well beyond the sensitivity of these superconducting devices. Our experiment in this case 
although unfeasible at present, is much more interesting from a theoretical point of view. Indeed for such dilute 
systems our basic assumption that the RE ions are independent from one another is completely justihed and therefore 
the theoretical interpretation of the results is much cleaner. Also experiments on gels, however, could become soon 
feasible improving the experimental layout e.g., with a smaller device-sample distance, improving the SQUID with 
a dedicated design, or in the next future using a free electron laser in order to have a photon flux in the ultraviolet 
energy range three to four orders of magnitude higher then that of a storage ring. 

We think that the successful detection of collective excitations modes such as scissors modes in RE compounds 
would trigger the search of a similar behavior in many other systems and, in particular, in actinides samples in which 
paramagnetic and deformed atoms can be present in different concentrations. 
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